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We discuss recent theoretical developments concerning the description of the 
production and decay of heavy quarks and colored scalars in the framework of 
nonrelativistic QCD. 
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1. Threshold Physics at the ILC 

The e + e~ center-of-mass (cm.) energy at a linear collider (LC) can be very 
precisely monitored, allowing for an accurate exploration of the threshold 
regime. The top-quark mass can be determined from a measurement of 
cr(e + e _ — > Z* 7 7* — » tt) line shape at a LC operating at cm. energies 
around the tt threshold (y/q* ~ 350 GeV). The rise of the cross section 
with increasing cm. energy is directly related to the mass of the top quark. 
Assuming a total integrated luminosity of 300 fb _1 , LC simulations of a 
threshold scan of the top-antitop total cross section have demonstrated 
that experimental uncertainties below 100 MeV for the top-quark mass 
determination can be obtained 1,2 , even when beam effects, which lead 
to some smearing of the effective cm. energy, are taken into account. If 
the normalization of the cross section line shape is well under control, it 
is possible to determine the strong coupling, the total top quark width 
and, if the Higgs boson is light, the top Yukawa coupling. In view of the 
accuracy obtainable at the LC the theoretical uncertainties for the total 
cross section should be lowered to a level of a few percent 3:4 . A precise 
knowledge of the top mass would also improve the analysis of electroweak 
precision observables and put indirect constraints on New Physics. It has 
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been shown 5 that an accuracy of 100 MeV on the top quark mass would 
allow to perform stringent internal consistency checks of the SM and of 
some scenarios of Supersymmetry (SUSY). 

Threshold studies would also be feasible for squarks at a next LC. Many 
models of SUSY predict that, due to large mixing, the lightest squark could 
correspond to one of the mass eigenstates of the third generation with 
nriq < 500 GeV, thus allowing for the production of stop pairs at a future 
e+e~ LC operating below 1 TeV. A "threshold scan" of the total cross 
section line-shape at such a facility will yield precise measurements of the 
stop mass, lifetime and couplings 6 , in close analogy to the program carried 
out in threshold studies for the top-antitop threshold 2 . 

2. Theoretical Status of it Production 

Close to threshold the top quark pairs are produced with small veloci- 
ties v <C 1 in the cm. frame. Therefore the relevant physical scales gov- 
erning the top-antitop dynamics are the top quark mass m t , the relative 
three-momentum p ~ nitv and the top quark nonrelativistic kinetic energy 
E ~ m t v 2 . Since the ratios of the three scales can arise in matrix elements, 
the cross section cannot be calculated using the standard QCD expansion 
in the strong coupling a s . The best known indication of the latter comes 
from the well-known "Coulomb singularity" , which shows up as a singular 
(a s /v) n behaviour in the v — > limit of the tt production amplitude at 
the n-loop order in perturbative QCD. The proper expansion scheme for 
the ti threshold region is a double expansion in both a s and v, and one 
has to use the parametric counting a s ~ v <c 1 to identify all effects con- 
tributing to a certain order of approximation. In this fixed-order expansion 
the leading order (LO) contributions correspond to terms in the total cross 
section proportional to v(a s /v) n , (n = 0, ...,oo), next-to-leading (NLO) 
terms are proportional to v(a s /v) n x [a s ,v] and so on. At LO the total 
cross section is proportional to the absorptive part of the Green function of 
a Schrodingcr equation containing the static QCD-potential 7 . Higher or- 
der corrections in this expansion are rigorously implemented employing the 
concept of nonrelativistic effective quantum field theories, first proposed 
by Caswell and Lepage 8 . In this scheme, the original QCD Lagrangian 
is reformulated in terms of an effective nonrelativistic Lagrangian called 
"Nonrelativistic Quantum Chromodynamics" (NRQCD) by using the hi- 
erarchy m t 3> p » E, which allows to separate short-distance physics at 
the "hard" scale of order the heavy quark mass from long-distance physics 
at the nonrelativistic scales p and E. The hard-momentum effects are en- 
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coded as Wilson coefficients of the operators in the effective Lagrangian. 
Operators with increasing dimension are introduced in the effective theory 
to include the effects of higher orders in the nonrelativistic expansion, but 
only a finite number is needed for a given precision. The NRQCD factor- 
ization properties were used in a number of NNLO calculations 9 ' 10 of the 
total ti production cross section. 

A common feature of the fixed-order calculations is that the running 
from the hard scale m down to the nonrelativistic scales was not taken into 
account. At NNLO NRQCD matrix elements and Wilson coefficients involve 
logarithms of ratios of the hard scale and the nonrelativistic scales, which 
in the case of top quark pair production close to threshold can be sizeable 
(for example a s (m t ) log[m t /E] ~ 0.8 for v ~ 0.15). Considering the para- 
metric counting a s log v ~ 1 introduces a modified expansion scheme for 
the size of the terms contributing to the total cross section, where the dom- 
inant contribution, proportional to (a s /v) n X)i( a s l°g u )% i n — 0, ...oo), 
is called leading- logarithmic order (LL). A number of different versions of 
NRQCD 11 > 12 ) each of which aiming (in principle) on applications in differ- 
ent physical situations, allow for rcnormalization group improved calcula- 
tions (see also Ref. 13 for recent reviews on the field). The EFT vNRQCD 
("velocity" NRQCD) 12 ' 14 has been designed for predictions at the ti thresh- 
old. It treats the case m t » p E > Aqcd, i-e. all physical scales are 
perturbative, but also has the correlation E t — pi/ m t built in at the field 
theoretical level. QCD effects for the ti total cross section up to the NNLL 
order computed in this framework 4 significantly reduce the size and scale 
dependence to yield a 6% theoretical uncertainty. 

The effective Lagrangian for the scalar version of vNRQCD, which de- 
scribes the nonrelativistic interaction between pairs of colored scalars and 
provides the needed ingredients for a summation of QCD effects at NLL 
order, has been given recently 15 . 

3. Finite Width and Electroweak Effects 

Up to now, no systematic and complete treatment of electroweak effects 
in the total cross section for top quark pair production has been achieved 
beyond the LO approximation. The large top width, being of the same 
order than the nonrelativistic energy, is essential in the description of the 
ti threshold dynamics. It was shown 7 that in the nonrelativistic limit the 
top-quark width can be consistently implemented by the replacement E — ► 
E + iT t in the results for the total cross section for stable top quarks. 
Although this replacement rule can accommodate some of the NLO and 
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NNLO electroweak corrections, a coherent treatment at the conceptual level 
requires the use of an extended NRQCD effective theory formalism. 

In the NRQCD/vNRQCD framework as long as one is not interested 
in any differential information of the top decay, finite lifetime corrections 
to the total cross section can be regarded as short-distance information to 
be encoded in the Wilson coefficients of the NRQCD Lagrangian and the 
NRQCD currents. As the particles involved in the electroweak corrections 
can be lighter than the top quark they can lead to nonzero imaginary parts 
in the matching conditions. These electroweak absorptive parts render the 
NRQCD Lagrangian non-hermitian, but the total cross section can still be 
obtained from the imaginary part of the e + e~ — > e + e~ forward scatter- 
ing amplitude by virtue of the optical theorem and the unitarity of the 
underlying theory, in complete analogy to the treatment of the inelastic 
processes in quantum mechanics where particle decay and absorption are 
implemented through potentials with complex coefficients. An important 
feature of the effective theory treatment of electroweak effects is that res- 
onant and non-resonant contributions in NRQCD amplitudes can be sys- 
tematically separated if the scaling relation v <~ a s ~ a 1 / 2 is used. The 
latter is justified because numerically the top width is approximately equal 
to the typical top kinetic energy T t ~ m t a ~ Ey^ n ~ m t a 2 s . Including the 
non-resonant background diagrams leading to the same final states as those 
of top decay in the matching calculations is necessary in order to maintain 
gauge invariance. 

In this approach the NNLL matching conditions accounting for the ab- 
sorptive parts related to the bW + bW~ final state were derived in 16 and 
shown to amount numerically to several percent. A very interesting new 
conceptual aspect of these corrections is that they have UV divergences 
that arise from the high energy behaviour of the tt effective theory phase 
space integration. The phase-space in the full theory is cut-off by the large 
top mass, but it extends to infinity in NRQCD where we have taken the 
limit rrit — > oo a . The divergences show up only when L t ^= 0, as a fi- 
nite width generates a distribution for the top invariant mass and thus 
allows for arbitrary large momenta in the nonrclativistic phase space in- 
tegration. Similar divergences had already been noted in the QCD NNLL 
relativistic corrections to the S'-wavc zero distance Green function if the 
unstable propagator was used 4 ' 9 as well as in the leading order P-wave 



a Instead of using the optical theorem to obtain the total cross section from the forward 
scattering amplitude one can integrate over the phase-space explicitly with cuts and 
these UV divergences obviously do not arise 
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zero distance Green function which accounts for tt production through the 
Z-exchange 4 . The NNLL divergences renormalize (e + e~)(e + e~) operators 
that contribute to the total cross section through the imaginary parts of 
their Wilson coefficients. The running induced in these Wilson coefficients 
by the divergences thus represents a NLL effect to the total cross section 16 
but their matching conditions at the hard scale are presently unknown. 

Scalar particles are produced at leading order in the nonrelativistic ex- 
pansion in a P-wave state. In these systems the phase-space divergences 
constitute a more severe problem as they show up already in the leading 
order Green's function. A phenomenologically well motivated example is 
pair production of SUSY partners of the quarks at threshold, which has 
been studied at LO in several works 17 using a semi-phenomenological so- 
lution in order to deal with the phase-space UV divergences. 

4. Nonrelativistic Currents with General Quantum 
Numbers 

There are a number of issues concerning the consistent formulation of non- 
relativistic interpolating currents which describe color singlet heavy quark- 
antiquark and squark-antisquark pair production for general quantum num- 
bers in n — 3 — 2e dimensions. These involve in particular the generalization 
of spherical harmonics in n spatial dimensions and the role of evanescent 
operator structures (i.e. that vanish as e — > 0) for the description of the 
nonrelativistic spin. The latter have been addressed in detail in a recent 
work 18 and we briefly comment on them here. 

Let us discuss first the spin singlet currents with arbitrary angular 
momentum L ( 2S+1 Lj — 1 Ll). The generic structure of the produc- 
tion currents with total spin zero is ^ p (x)T(p)x*_ p (x), where T(p) rep- 
resents an arbitrary tensor depending on the cm. momentum label p and 
X*L p = (i^2)X*~ P - The interpolating currents associated to a definite angular 
momentum state L are related to irreducible representations of the tensor 
F with respect to the rotation group SO(n). The irreducible tensors are up 
to normalization just the spherical harmonics Ylm(ti,p) of degree L, with 
M = 1, . . . ,til, that form an orthogonal basis of a nL-dimensional space 
with Ul = (2L + n — 2) r(^-i)r("£+i) • ^ representation in terms of cartesian 
coordinates of the spherical harmonics of degree L is given by the totally 
symmetric and traceless tensors with L indices T tl — %L (p), where the indices 
i\ . . . , iL are cartesian coordinates: 

T n-iL {p) =p h^ p i L _ P (gfafa i, .. . p <r. + ...)+.. . , (1) 

IL + n — 4 
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and which satisfy the eigenvalue equation for the squared angular momen- 
tum operator L 2 T h - iL (p) = L(L + n - 2) T il - iL {jp). The currents with 
angular momenta S, P and D are e.g. relevant in the electromagnetic pro- 
duction of colored scalars from e + e~ and 77 collisions. The use of the 
generalized currents built from (1) is mandatory to obtain consistent re- 
sults in dimensional rcgularization in accordance with SO(n) rotational 
invariance. An instructive example regarding the computation of the non- 
relativistic three-loop vacuum polarization diagram with two insertions of 
the Coulomb potential is explained in Ref. 18 . 

The interpolating currents describing the production of a fermion- 
antifermion pair in a spin triplet S = 1 state for arbitrary L ( 3 Lj) re- 
quires the treatment of Pauli cr-matrices in n-dimensions. The er-matrices 
a 1 (i = 1, . . . ,n) arc the generators of SO(n) rotations for spin 1/2 and 
satisfy the Euclidean Clifford algebra {a 1 , cr-?} — 25 %3 . As for the case of the 
7-matrices 19 products of cr-matrices in arbitrary number of dimensions 
cannot be reduced to a finite basis, but represent an infinite set of inde- 
pendent structures, which can be choosen as the antisymmetrized product 
of <7 matrices: a h - im = a^a* 2 ■ ■ ■ a 1 ^ , m = 0, 1, 2, ... For m < 3 the 
<7 n "'* m are related to physical spin operators, with eigenvalues of S 2 equal 
to (0, n — 1, 2, n — 3), which reduce to the known n = 3 values. The to > 3 
operators are evanescent for n ^ 3 (although their spin eigenvalues are 
non-zero). S-wave currents in an arbitrary spin state thus have the form 
ipp{x) a 11 '" lm X*- P ( x )> an d can arise in important processes. The structure 
a tlt2t3 for example arises in fermion pair production in 77 collisions, while 
the evanescent operator a 11 "" 15 is present in the // — > 37 annihilation ampli- 
tude. Note that the differences between the two different singlet (m = 0, 3) 
and triplet (to = 1,2) currents correspond to evanescent operators as well. 

It is well known from sublcading order computations based on the ef- 
fective weak Hamiltonian that one needs to consistently account for the 
evanescent operator structures that arise in matrix elements of physical 
operators when being dressed with gluons. A renormalization scheme can 
be adopted such that a mixing of evanescent operators into physical ones 
does not arise 19 . Moreover it is also known 20 that modifications of the 
evanescent operator basis correspond to a change of the renormalization 
scheme. While this does not affect physical predictions, it does affect ma- 
trix elements, matching conditions and anomalous dimensions at nontrivial 
sublcading order. Thus precise definitions of the schemes being used have 
to be given to render such intermediate results useful. 

In the framework of the nonrelativistic EFT these properties still apply. 
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However, using the velocity power counting in the EFT allows for even more 
specific statements. Concerning interactions through potentials, transitions 
between the different S-wave currents built from a ll '" %m cannot occur be- 
cause the potentials are SO(n) scalars and the currents are inequivalent 
irreducible representations of SO(n). Even for currents with L ^ and for 
the spin-dependent spin-orbit and tensor potentials (which is all we need 
to consider at NNLL order) one can show that transitions between cur- 
rents containing cr ll '" %m with a different number of indices cannot occur 18 . 
The same arguments apply to the exchange of soft gluons in vNRQCD. 
Concerning the exchange of ultrasoft gluons, transitions between currents 
containing cr* 1 '"*"* with a different number of indices can arise, but only if 
the interaction is spin-dependent. The dominant among these interactions 
corresponds to the operator tp^a 13 k 3 ip p A l and can only contribute at N 4 LL 
order, which is beyond the present need and technical capabilities. 

So for the S-wave currents in n dimensions one can employ either one 
of the two spin singlet (m = 0, 3) or triplet currents (m = 1,2) in the EFT 
and the difference corresponds to a change in the renormalization scheme. 
This means in particular that as long as the renormalization process is 
restricted to time-ordered products of the currents, one can freely use three- 
dimensional relations to reduce the basis of the physical currents. However, 
once the basis of the physical currents is fixed, one has to consistently apply 
the computational rules in n dimensions. Moreover, one can also conclude 
that currents containing evanescent o %1 '" %m matrices (m > 3) can be safely 
dropped from the beginning as long as one does not need to account for 
spin-dependent ultrasoft gluon interactions. 

Based on the latter considerations it is straightforward to construct 
spin-triplet currents with arbitrary L ( 3 Lj). They can be obtained 18 by 
determining irreducible SO(n) representations from products of the tensors 
T %1 '" %L (p) describing angular momentum L and the spin-triplet S = 1 
currents discussed previously. As for the case of the S-wave currents the 
physical basis for arbitrary spatial angular momentum is not unique due to 
the existence of evanescent operator structures. It is possible 18 to construct 
currents with fully symmetric indices equal to the total angular momentum 
J by using the two spin triplet operators, a 1 and a %3 . There are also currents 
having more than J indices, which transform according to more complicated 
patterns of SO(n). They become equivalent in n = 3 to the fully symmetric 
currents 18 and are also appropriate to describe production of 3 Lj states. 

The NLL anomalous dimensions of the currents with arbitrary spin and 
angular momentum configurations have been calculated in Ref. 18 . Since 
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at LL order the currents are not renormalized, their NLL order anomalous 
dimensions are independent of the scheme used for the currents or the 
potentials. The NLL running found for the currents shows a suppression 
oc 1/(2.L + 1), which suggests that the summation of logarithms of v for 
the production and annihilation rates of high angular momentum states is 
less significant. 
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